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ABSTRACT 

The evolution of baryon density perturbations on very small scales is investigated. 
In particular, the nonlinear growth induced by the radiation drag force from the shear 
velocity field on larger scales during the recombination epoch, which is originally 
proposed by Shaviv in 1998, is studied in detail. It is found that inclusion of the 
diffusion term which Shaviv neglected in his analysis results in rather mild growth 
whose growth rate is <C 100 instead of enormous amplification ~ 10 4 of Shaviv's 
original claim since the diffusion suppresses the growth. The growth factor strongly 
depends on the amplitude of the large scale velocity field. 

The nonlinear growth mechanism is applied to density perturbations of general 
adiabatic cold dark matter (CDM) models. In these models, it has been found in 
the previous works that the baryon density perturbations are not completely erased 
by diffusion damping if there exists gravitational potential of CDM. With employing 
the perturbed rate equation which is derived in this paper, the nonlinear evolution 
of baryon density perturbations is investigated. It is found that: (1) The nonlinear 
growth is larger for smaller scales. This mechanism only affects the perturbations 
whose scales are smaller than ~ 1O 2 M0, which are coincident with the stellar scales. 
(2) The maximum growth factors of baryon density fluctuations for various COBE 
normalized CDM models are typically less than factor 10 for 3 — a large scale velocity 
peaks. (3) The growth factor depends on fib- 

Subject headings: cosmology: theory 

1. INTRODUCTION 

Structure formation of the universe from the large scale structure to the stellar objects is 
one of the most important issues of modern cosmology and astrophysics. Within a framework of 
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the hierarchical clustering scenario which is suggested by the adiabatic cold dark matter (CDM) 
models, formation of first luminous objects in the universe has been discussed by many authors 
(e.g., Fukugita & Kawasaki 1994; Pnedin fc Ostriker 1997 ; Ostriker &: Gnedin 1996| ; Haiman &; 



Loeb 1997; Nishi, et al 1998). According to those theoretical investigations, in which the initial 
dark matter and baryon density fluctuations have been solved linearly, the first luminous objects 
are thought to appear in the gas clouds in the dark halo at z ~ a few x 10, depending on the 
cosmological models. 

On the other hand a new mechanism which causes rapid growth of the baryon density 
fluctuations on very small scales during a short period of the recombination epoch has been 
proposed by Shaviv (1998). This growth of the baryon density perturbations is caused by the 
radiation drag force from the large scale velocity field which is coupled to the small scale baryon 
perturbations in a nonlinear manner. He has pointed out enormous growth (~ 10 ) of baryon 
fluctuations if the wavelength of the perturbations is shorter than a critical value. Such enormous 
amplification of baryon density fluctuations must lead to early formation of stellar objects. If the 
amplification is as large as Shaviv's prediction, therefore, it might be a quite interesting antithesis 
for the standard scenario of the cosmic structure formation in the CDM cosmological models, 
though we will show that this nonlinear growth mechanism would not be so drastic. 

In this paper, we first carefully re-examine Shaviv's treatment of density perturbations. 
We find the diffusion term which Shaviv has neglected in his analysis plays a crucial role for 
suppression of the nonlinear growth. 

There are also some artificial assumptions in Shaviv's investigation. First one is in use of 
the thermal equilibrium for the recombination process. This nonlinear growth is induced by the 
shear velocity field on larger scales through the radiation drag force against density fluctuations of 
ionized hydrogens. Therefore we need to calculate fluctuations of the ionization fraction. Shaviv 
has obtained these fluctuations by employing the perturbed Saha equation which should not be 
applied during the recombination epoch. Here we develop a new method to treat these fluctuations 
by the perturbed treatment of the rate equation which is an extension of the previous work by 
Yamamoto, Sugiyama, and Sato (1997, Paper I). 

The next assumption Shaviv has made in his investigation is the existence of isothermal 
baryon density fluctuations. It is well known that the adiabatic baryon and photon density 
perturbations are erased away due to the photon diffusion before and during the recombination 
epoch (Silk 1968; Sato 1971; Weinberg 1971). The diffusion (Silk) damping scale goes over the 
galaxy scale. Therefore it seems to be unavoidable for the nonlinear growth mechanism to assume 
small scale isothermal baryon perturbations during the recombination epoch as the seeds of 
baryon perturbations in an artificial manner. The origin of the small scale isothermal baryon 
perturbations has not been considered. 



However we would like to stress that the small scale isothermal baryon density fluctuations 
are produced after the Silk damping in general CDM cosmological models without any ad hoc 
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assumptions. Such isothermal baryon density fluctuations on small scales have been found by 
Yamamoto, Sugiyama & Sato (1998, Paper II). Once baryon fluctuations are damped by diffusion. 
However they grow again owing to the gravitational potential of CDM before the recombination 
epoch due to the breakdown of the tight-coupling between baryons and background photons. Since 
the Compton drag from background photons is still effective at that epoch and works as friction, 
this growth of baryon density fluctuations is characterized by the terminal velocity. Such natural 
existence of the small scale isothermal baryon density fluctuations motivates us to revisit Shaviv's 
investigation in the context of the conventional adiabatic CDM models. 

The paper is organized as follows: In §2, we review the nonlinear growth mechanism which 
Shaviv has proposed. Some numerical tests for a toy model with simple assumptions are presented 
in §3. It is shown that the diffusion term plays a crucial role for the estimate of the growth 
factor. In §4, we examine whether the nonlinear growth mechanism is effective on the adiabatic 
CDM models or not. §5 is devoted to summary. Throughout this paper we work in units where 
c = fi = ks = 1, and assume To = 2.726 K as the cosmic microwave radiation temperature at 
present. 



2. Basic Equations for the Nonlinear Effect 

Before describing the nonlinear growth mechanism of small scale baryon density fluctuations 
proposed by Shaviv (1998), let us first summarize basic equations (see also Paper I; II). We work 
in the Newtonian perturbed metric, 

ds 2 = (^ 2 (-(l + 2^)d?? 2 + (l-2^)dx 2 ), (1) 

where ^ is the perturbed gravitational potential, a is the scale factor whose suffix indicates 
the present value, and rj = J dt(ao/a). As we are interested in the small scale cosmological 
perturbations, we can assume the geometry of the universe to be flat. Then, in the high redshift 
universe, the Friedmann equation leads to 

H 2 =( a -^) 2 =(^X^±^n Hl (2) 
\a a ) \a J a cq + a 

where the dot denotes ^-differentiation, Hq is the Hubble constant, and a cq is the scale factor at 
the matter-radiation equality epoch which can be written as 

^- = 4.04 x 10 4 (1 - f u )n h 2 , (3) 

with /„ being the neutrino fraction in the energy density of radiation components. For the 
massless standard neutrino model with three species, f u ~ 0.405, which we adopt throughout the 
present paper. Here, we take into account the fraction of helium abundance. Hence we define the 
fractional ionization 

*e = 7^ , (4) 

2n He + n H 
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where n e , n^ e , and nn are the number densities of free electron, helium, and hydrogen, respectively. 
The primordial helium mass fraction is defined as 

4n He 

V P = , (5) 

n b 

where n b = 4nn e + is the net baryon number density. We set y p = 0.24 throughout 
the present paper. The fractional ionization x e is 1 in the very early universe and becomes 
(1 — y p ) /(l — Hp/2) ~ 0.86 just before the recombination epoch due to the early recombination of 
helium atoms. 

Let us now describe the nonlinear growth mechanism. We proceed with reviewing equations 
of the small perturbations for the baryon-electron fluid. Since we are interested in the evolution of 
density perturbations on very small scales, we adopt the Newtonian approximation. The evolution 
equations for the baryon-electron fluid coupled with the primeval radiation through Compton 
scattering are written as 

Pb + 3- Ph + V i (p h V£) = 0, (6) 
a 

K + -K + + -v*p + v% = -^(v; - v£), (7) 

where p b an d are the density and the velocity of the baryon-electron fluid, respectively, P is 
the pressure, ot is the Thomson cross section, VI is the dipole moment of the photon field, and R 
is defined as R = Sp^/Ap^ with the spatially averaged energy density of the photon field p 7 , where 
~ denotes the background quantity. The right hand side of equation ([?]) represents the radiation 
drag force caused through the Compton interaction between electrons and photons which can be 
written by using td that is the time scale of momentum transfer between the baryon-electron fluid 
and the photon fluid as 



, ., , ./ 7 -K b . = — (vj-K), (8) 
ao orrip \ 2 J V ' / td 

where m p denotes the proton mass. As the universe expands, the momentum transfer time scale 
td becomes longer. However the drag time scale is still small enough to keep the tight coupling 
between baryon-electron and photon fluids even just before the recombination epoch, so that 
V* = V^. Eventually the tight coupling breaks down as the recombination process proceeds when 
the photon mean free path becomes larger than the wavelength of the density perturbations (see 
also Paper II). 

Let us now focus on perturbations of the baryon-electron fluid on very small scales. As far as 
considering linear density perturbations, perturbations of each scale evolve individually. Following 
Shaviv (1998), however, we take into account the quasi-nonlinear contribution of the drag force 
from the large scale velocity field. To be specific, we separate the velocity field into the small scale 
part which we focus on and the large scale part which describes the nonlinear contribution to the 
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small scale perturbations as, 

v; - v£ = v; (s) - k (s) + v; (L) - v^ (L) , (9) 

where indices (S) and (L) indicate small and large scales, respectively. The fluctuation of the 
fractional ionization on the small scale is defined as Sx e rg\ = x e — x e . Then the right hand side of 
equation (^) leads to 



(\ - 1) (x e (y; {s) - K (s) ) + fee(s)A^), (10) 



ao Snip 
where 

^^^(L)-^(L)- (11) 

Here we ignore higher order terms of the small scale perturbations. The term <5x e (g)AVQ in 
equation ([k]) describes the quasi-nonlinear coupling between the large scale shear velocity and the 
small scale baryon density fluctuations. Adding this quasi-nonlinear term, we can write the linear 
perturbation equations on the small scale as 

5b(S)+ViVj (s) =0, (12) 
K ( s) + ^( S ) + ^A(s)+Vit( S) = ±?^(-x e Vfa +5x e{s) Avi), (13) 

where c s is the sound velocity of the baryon-electron fluid defined by c s 2 = P/pb- In deriving 
equation (|l3|) , we have assumed the isothermal perturbations for the small scale perturbations, 
i.e., V\g\ = since the Silk damping process erases the density perturbation and the dipole 
moment of the photon field on small scales. 

Let us investigate the effect of the quasi-nonlinear term fe e (g)AVg on the evolution of 
perturbations. Before the recombination epoch, AVq = due to the tight coupling. However, the 
tight coupling becomes broken down, i.e., AVq / 0, as the decoupling process proceeds and the 
contribution from the quasi-nonlinear term may not be negligible anymore. In fact Shaviv has 
claimed that this quasi-nonlinear contribution, which represents the spatially fluctuating radiation 
drag force, causes enormous amplification of baryon density fluctuations during the recombination 
epoch. Hereafter, we omit the label '(£>)', for simplicity. Combining equations (^) and (|l3|), we 
obtain 

4 + ^ b - c s 2 A5 b - A* = -±5 h (l + Ay ° V / Xe ), (14) 
where A is the Laplacian. 

By assuming the wave form solution, 5^ oc e ( w? ?- 4x k ) ) where k is the comoving wave number, 
we have the dispersion relation 

- 2 +(4^M4-?- : 4H' (15> 
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where 



— = - = 3.2 x lO^V^o/iVl + z s _1 , (16) 
m a 

J_ = ^AOTP^ ( l _y JL \ =7Ax 1Q -24- e(1 + zf B -l (17) 

r D a 3m p V 2 / 

1 / Af \ 

_ = Cs k = 1.61 x l(T 16 (ft b /i 2 ) 1/3 ti- + + 2 s -1 , (18) 

\M Q J 

1 / ~~/ a \2 /30 b 1 . , 

rj V Va o y V 2^o t h 



J ^7l«l A ^fc_ 49xlo -17 (ie | a | Ayo) l/2 ( ^ 2) l/6 (1 + z) 3/2^\ 1/6 g -l 5 (2Q) 



~!( \ TD ' V M G 

with k = |k| and AVo = AVq • k l /k being the shear velocity from the large scale. Here th, Td> 
t , rj, and tr are time scales for the Hubble expansion, the radiation drag, the sound oscillation, 
the Jeans oscillation, and the quasi-nonlinear radiation force, respectively. The parameter a of 
equation ( p0[ ) is defined as 

(21) 

x e d h 

A nonzero value of a provides the quasi-nonlinear effect. Shaviv obtained a by employing the 
Saha equation. The validity of this assumption will discuss in §4. The derivation of 5x e is shown 
in Appendix B. In the above expressions, the baryon mass M is used instead of the comoving 
wave number k (Paper I), which is defined by 

M = j*P*f«*\\ (22) 
3 \kaj y 1 

It should be noticed that we take into account only the baryon component for the gravitational 
potential. 

Let us now evaluate the growth rate of the growing-mode solution. The solution of the 
dispersion relation (|l5|) can be written as 



1 

W= 2 



i i y j i i .i 



1 1 

Th Td/ V V t h T D/ \ T o r J r R 



(23) 



Here we omit the negative solution which provides the decaying mode solution. The positive 
real part of u, i.e, K[o;(ry)], represents the growth rate of the instability. The exact expression of 
K[cj(r7)] is shown in Appendix A. In the time interval between 771 and 772, the amplification of the 
baryon density fluctuations, which we refer the growth factor D (772,7/1), becomes 

D( V 2,m) = ^ = exp( £ d V M[u( V )\) . (24) 
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3. Nonlinear Growth of a Toy Model 

Before carrying out fully consistent calculations based on density perturbations of actual 
cosmological models, we investigate the nonlinear growth for a toy model. First, we assume the 
large scale velocity AVo is constant in time throughout this section. Secondly, as Shaviv has done 
in his paper, the fractional ionization x e and its perturbation 5x e are calculated by employing 
the Saha equation, which is described in Appendix B. We take the fiducial CDM model, i.e., 
= 1.0, fib = 0.04, and h = 0.5. The time evolution of the fractional ionization is shown in 
Fig. 1. We also plot the redshift visibility function 

«,) - (25) 

where r = — Jjj° drf — n e OT is the optical depth. The width of the visibility function corresponds to 
the recombination epoch. Therefore the recombination epoch of this model is 1 + z = 1500 — 1100. 

Let us first compare time scales of equations (|T^)-(p0|). In Fig. 2, these (inverse) time 
scales are plotted as a function of (1 + z). The panel (a) shows 1/td, 1/th, and 1/rj, all which 
do not depend on k. Note that 1/td dominates others before and during the recombination 
epoch. However as the recombination process proceeds, this time scale suddenly increases and 
eventually exceeds other time scales after recombination. In this panel, we plot t /t^ which is 
also independent on k. We will discuss this variable later in this section. 

In the panel (b), we plot l/r Q and 1/tr for different wave numbers k, i.e., k = 71.2, 330, 
1530,7120, and 33000MPC" 1 which correspond to mass scales 10 6 , 10 4 , 10 2 , 10°, and 1O~ 2 M , 
respectively. 

Now let us discuss the growth rate 3?[a;(r/)] by comparing the time scales in the expression 
(|23|). As is shown in Fig. 2, the larger the wave number k, the larger the inverse time scales 1/tr, 
and 1/t become. Therefore we expect to obtain the maximum growth rate when we take the 
small scale limit, i.e., k — > oo. As is shown in Appendix A, the real part of the square root term of 
equation ( |23| ) approaches to t q /t^ which is the combination seen in the panel (a) of Fig. 2. Shaviv 
claimed that this combination provides the maximum growth rate, i.e., 



However, this is not the whole story. The complete expression of equation ( [23] ) with large k 
limit becomes 



RMiflW = Urn = --- + ^ = --- + ^ . (27) 



1J_ To_ _ 1 1 | aAV 
k— >oo 1 ' 2 tq 2t^ 2 td 2c s td 

Here we ignore the inverse Hubble time scale because it is much smaller than 1/td during 
recombination. Remember that 1/td is comparable to t q /t^ during recombination as is shown in 
Fig. 2(a). In fact, the existence of this 1/td term suppresses the nonlinear growth a lot. From 
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Fig. 1. — Fractional ionization x e (dashed line) and visibility function g(z) (solid line) from the 
Saha equation with the CDM model; f2o = 1.0, $\ = 0.04, and h = 0.5. The peak of the visibility 
function (z = 1300) corresponds to the recombination epoch. In this fiducial CDM model, the 
recombination occurs in between 1 + z = 1500 and 1100. 
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Fig. 2. — Inverse time scales (eq. [p| - eq. ^0[) for the fiducial CDM model as in Fig. 1. The 
Saha equation is employed to calculate the fractional ionization and AVo = 10~ 4 for tr. Wave 
number independent and dependent time scales are plotted in (a) and (b), respectively. In panel 



(a), t q /t^ which provides the maximum growth rate together with td (eq. |27|) is also plotted by a 
long dashed-short dashed line. In panel (b), the wave numbers of the baryon density perturbations 
are k = 71.2,330,1530,7120, and 33000MPC" 1 which correspond to mass scales 10 6 , 10 4 , 10 2 , 10°, 
and 1O _2 M0, respectively, from bottom to top. 
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equation (|27]), it is found that a sufficient condition for the nonlinear growth at least on the 
smallest scale is K[w(7y)] max > 0, i.e., 

aAV > c s . (28) 

In Fig. 3, the growth rates and the growth factors of the toy model are plotted for various 
wave numbers and the large wave number limit. 

From this figure, it is found that the growth rate and the growth factor are very sensitive to 
the value of the large scale velocity field AVo- The growth factor of AVq = 10 -4 is more than 
factor 20 larger than the one of AVo = 6 x 10~ 5 for the large k limit. 

The "maximum" growth rates and growth factors obtained by Shaviv's formula (eq. p6[) 
are also plotted for comparison in Fig. 3. It is clear that neglect of the diffusion results in the 
extremely large false growth factor. By the correct treatment, rather mild growth is found. When 
we take relatively large velocity AVb = 10 -4 , we obtain about factor 50 growth factor while 
Shaviv's treatment induces more than 10 5 growth of small scale density perturbations. From 
Fig. 3, it is shown that the nonlinear growth mechanism starts to work in the relatively early 
epoch (z ~ 1500) if we neglect the diffusion term. However, by the correct treatment, it is found 
that the nonlinear growth works only when t g /t^ > 1/td, i.e., 1 + z < 1340 in this fiducial CDM 
model as is shown in Fig. 2. Even in this shorter period, the maximum growth rate is suppressed 
by the existence of the diffusion term (eq. [^7]]). We should conclude that the enormous growth 
~ 10 4 of the baryon density fluctuations which Shaviv claimed in his paper is not likely to happen. 

However, we should emphasize that this nonlinear growth mechanism still works and enhances 
the baryon density perturbations. It is interesting to estimate the nonlinear growth factor for 
realistic cosmological models. 

4. Nonlinear Growth in adiabatic CDM Models 

Since it is generally believed that adiabatic small scale baryon density perturbations decay 
out before and during the recombination epoch due to the diffusion (Silk) damping, an additional 
assumption to standard adiabatic perturbations seems to be needed for the nonlinear growth 
mechanism which requires the existence of small scale baryon perturbations. In the paper by 
Shaviv (1998), for example, he artificially assumed the existence of baryon isothermal density 
perturbations. 

However, we should point out that small scale isothermal baryon density perturbations are in 
fact existing at the recombination epoch in general adiabatic CDM models. Recently, the authors 
of the present paper (K.Y. & N.S.) and Sato have found that small scale baryon perturbations 
grow after the diffusion (Silk) damping even before the recombination epoch (Paper II). This 
growth of baryon density perturbations before the recombination epoch is the result of the 
breaking of the tight-coupling between the baryon perturbations and the photon perturbations 
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Fig. 3. — Growth rates ^R-[uj(r](z))] ((a) and (c)) in the unit of the Hubble time and growth 
factors D(t](z), 0) = exp (fg drf^t[u{rf)]) = exp (/ 1 °^ 2 d(l + z')M[u]/H^ ((b) and (d)) of two simple 
simulations with fixed shear velocities AVo = 10 -4 and 6x 10 -5 for the same model as in Fig. 1. The 
dotted lines plot the growth rates or growth factors for the wave numbers k = 33000, 7120, 1530, 330, 
and 71.2Mpc _1 from top to bottom. One can see almost no growth for k = 330 and 71.2Mpc . 
The solid lines are the maximum growth rates (eq. pTfl ) and growth factors. The maximum growth 
rates and growth factors of the investigation by Shaviv (eq. (2(| ) are also plotted by dashed lines for 
comparison. It is shown that inclusion of the diffusion term 1/td (eq. |27j) results in the significant 
suppression of the growth factors. 
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on small scales. The growth is characterized by the terminal velocity, which implies the balance 
between the gravitational force due to the potential of CDM and the friction force by the Compton 
interaction with background photons. According to the previous investigation in Paper II, the 
Fourier coefficient of the small scale baryon perturbations grows in proportion to (1 + z)~ 7//2 . The 
small scale baryon density perturbations have amplitude of order O(l0~ 2 ) ~ O(10~ 3 ) compared 
with the CDM perturbations. 

The small scale photon perturbations are smeared out by the diffusion damping. On the 
other hand, the energy transfer between the baryon-electron fluid and the background photons 
is effective during and for some time after the recombination epoch. Therefore the nature of the 
small scale baryon perturbations is isothermal during these epochs (see Paper I) . These isothermal 
baryon perturbations could be enhanced by the nonlinear growth mechanism although the growth 
factor may not be as large as originally thought by Shaviv as is described in the previous section. 

The amplitude of the large scale velocity field AVo is crucial for the instability as is shown in 
Fig. 3. In linear theory of density perturbations, the variance of square of the large scale velocity 
field is evaluated by 

(1^ _ y b f ) = _L Jdkk 2 \v^(k, v y-v h (k, v y\ 2 , (29) 

where V^ik^rf) 1 and Vb(/c,r/) J are Fourier coefficients of the dipole anisotropy of the photon field 
and the peculiar velocity field of the baryon fluid, respectively. Introducing the parameter u, we 
derive the large scale velocity field from the rms one as 

AV = ■ (30) 

We consistently solve the linear perturbations of the photon, baryon, and dark matter system 
in the expanding universe to calculate the rms velocity on large scales. In Fig. 4, the rms velocities 
of various COBE normalized CDM models are shown. It is found that the rms velocities are 
typically 10 -6 just before recombination. They rapidly increase toward 10 -4 as the recombination 
process proceeds when the tight coupling between photons and baryons breaks. We may conclude 
that effective values of the rms velocities during the recombination epoch are about a few times 
1CP 5 regardless of model parameters. According to our toy model calculations as is shown in 
Fig. 3, therefore, we expect only high velocity peaks induce the nonlinear growth on small scales 
in the realistic cosmological models. 

To calculate the growth rate and the growth factor of baryon density fluctuations, we need to 
obtain perturbations of the fractional ionization 8x e . Once we can get 5x e , we calculate a, which 
is defined by equation ([H]) to estimate tr. Using baryon density fluctuations 5b of the previous 
time step, we can calculate tr, the growth rate and the growth factor. We can estimate <5b with 
multiplying this growth factor by the amplitude of linear density perturbations with the COBE 
normalization. 



-13- 




2000 1500 1000 500 

1+z 



Fig. 4. — Root mean square shear velocities in large scales. The COBE normalized velocities of 
various flat CDM models grow typically from 10~ 6 to 7 x 10~ 5 during the recombination epoch. 
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Strictly speaking, we may have to solve radiative transfer during the recombination epoch to 
get precise value of fluctuations of 8x e . Instead, we have employed two approximation methods 
here. 

First method is to employ Saha's formula as Shaviv has done in his paper, which is a very 
simple but unrealistic assumption during and after the recombination epoch. It is well known 
that even the fractional ionization itself leaves from the value by Saha's formula as recombination 
proceeds (see e.g., Jones & Wise 1985). Therefore we next develop a perturbed treatment of the 
rate equation. Detailed treatments of these two methods are described in Appendix B. It is shown 
there that these two methods provide very different x e and 5x e . 

It is found that the nonlinear growth mechanism does not work at all for the fiducial CDM 
model with = 1, S\ = 0.04, and h = 0.5 if we employ the Saha equation for the evolution of 
x e and 5x e . It is because AVo is too small during the recombination epoch, i.e., 1 + z = 1500 and 
1100. In this calculation, 1/th + 1/td > t q /t^ for all time, which leads to the real part of uj of 



equation (23) to be negative on small scales. 



If we employ the perturbed rate equation to calculate Sx e , however, the nonlinear growth 
mechanism is effective since the recombination process continues until much later epoch, z ~ 800 
as is shown in Fig. 7 in Appendix B, when AVo is large enough to induce the nonlinear growth 
by the drag force. We obtain the nonlinear growth factor as a function of redshift z in Fig. 5 for 
various COBE normalized CDM models, i.e., the fiducial model with Oo = 1 and h = 0.5 and low 
density flat and open models with S7o = 0.3 and h = 0.7. For the purpose of comparison, different 
values of baryon density, i.e., Q^h 2 = 0.01 and Q^ii 2 = 0.02 are considered. We adopt 3 — a peaks 
for the amplitude of the large scale velocity field. The nonlinear growth mechanism works from 
z = 1100 to z = 800 ~ 600 which depends on the value of fl^h 2 . 

Fig. 6 plots the final growth factors as a function of wave number k. We find that the growth 
factors of low density CDM models are less than the ones oi Qo = 1 models. It is mostly because 
the COBE normalized rms velocity fields of low density models are smaller when the mechanism 
effectively works (z = 1100 ~ 800) as is shown in Fig. 4. The strong dependence of the growth 
factor is also found in Fig. 6. The reason is following. The residual ionization at low redshifts 

1 /2 

scales clS (X Q ' /(Q h h) (see e.g., Peebles 1993). Therefore we have smaller x e for larger Qb- The 



maximum growth rate can be written as (eq. |27|]), K[o;(^)] max = (aAVo/2c s — 1) / (2td) . As is 



shown in Fig. 8 in Appendix B, the fib dependence of a is relatively weak because of cancellation 
in between x e and 5x e . The remaining dependence appears through 1/td oc x e . For larger f?b, 
therefore, x e and the maximum growth factor become smaller. 



5. Summary 



We investigate the evolution of the baryon density perturbations on very small scales, 
especially focus on the nonlinear growth induced by the shear velocity field on large scales, which 
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Fig. 5.— Growth factor D(rj(z),0) = exp^JJ dr/9t[w(rf)]J for various CDM models. The large scale 
velocity filed AVo is consistently solved by linear density perturbations with COBE normalization 
(see Fig. 4). The 3 — a velocity peaks {y = 3) are considered. The fractional ionization and its 
fluctuation are computed by the rate equation. In parentheses, the cosmological parameters h, Qq, 
and f^b are shown. The mass scales of the baryon perturbations are 10~ 2 , 10°, 10 2 , 10 4 , and 
1O 6 M from top to bottom, respectively. Panel (b) and (e) are low density flat models while (c) 
and (f) are open models. 
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Fig. 6. — Final growth factor D(rj(z = 100), 0) for various cosmological models with v = 3 as a 
function of the wave number. The growth factor of the smallest baryon perturbation scale we chose 
(M = 1O~ 2 M ) is close to the maximum growth factor exp(j dr/K[u;(r/)] max ^ , where 5ft[w(r/)] max is 
defined in equation (27). 
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is originally proposed by Shaviv (1998). 

First we study this nonlinear evolution with employing simple assumptions, i.e., the constant 
large scale shear velocity field and the Saha equation for the calculations of the fractional 
ionization and its fluctuation. What we find is not enormous amplification of density fluctuations 
but relatively mild growth. The reason of Shaviv's overestimation is because he has neglected the 
diffusion term which suppresses the growth. 

It is found that the growth factor is very sensitive to the value of the large scale shear velocity 
field. Only small difference (AVb = 10 -4 versus 6 x 1CP 5 ) causes huge divergence of the final 
growth factor (factor 50 versus 2.5). 

Following the previous works (Paper I &: II), we next apply the nonlinear growth mechanism 
to density perturbations of general adiabatic CDM models. In these works, it has been found that 
the baryon density perturbations are not completely erased by diffusion damping if there exists 
gravitational potential of CDM. They even grow before recombination under the balance between 
the radiative drag force and the gravitational force. We investigate the nonlinear growth of these 
baryon density fluctuations with employing the perturbed rate equation which is proposed first 
time in this paper. The results are followings: (1) The nonlinear growth is larger in smaller scales. 
This mechanism only affects the perturbations whose scales are smaller than 10 2 Mq. (2) The 
maximum growth factors of baryon density fluctuations by this mechanism for various COBE 
normalized CDM models are typically less than factor 10 even if we take 3 — a large scale velocity 
peaks. (3) The strong dependence of the growth factor is found. This is because of the 
dependence of the residual ionization. The fractional ionization x e , which is proportional to the 
inverse diffusion time scale and the maximum growth rate, is larger for smaller f^- 

How could this nonlinear growth mechanism affect on the structure formation in the high 
redshift universe? It is very interesting coincidence that the typical scale of the nonlinear growth 
is about the stellar size 1O 2 M ). However we need to carry out high resolution numerical 
simulations to investigate the evolution of the stellar scale density fluctuations on which a lot of 
other complicated physical processes, i.e., shock heating, UV radiation from first starts or AGN, 
cooling and so on, work. Our quasi-linear perturbation analysis provides the initial condition 
for the calculations of such complicated structure formation in the high-z universe. The transfer 
function of the baryon density perturbations can be obtained by multiplying the nonlinear growth 
factor by the transfer function of the linear perturbations which is derived in Paper II. 
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A. Real Part of the Complex Square Root 

The derivation of the positive real part of the complex square root term in equation ( |T5| ) 
is shown in this appendix. In general, the square root of an arbitrary complex number can be 
rewritten as 



V 'A + Bi = \] V A 2 + B 2 (cosa + i sin a) = \j \J A 2 + B 2 exp ( y ) , (Al) 



where A and B are real numbers and 



A . B 

cos a = = , sin a = =. (A2) 

VA 2 + B 2 VA 2 + B 2 V ' 



Hence the real part of this square root is 



U[VA + Bi] = \JVA 2 + B 2 cosa/2 



Let us now take 



Ua+^A 2 + B 2 ). (A3) 



/ 1 1 \ 2 1 1 

A = - + - + 4 -" 4 -' ( A4 ) 
\Th to J rj 

5 = 4^-. (A5) 



R 



In the small scale limit, i.e., k — ► oo, we have l/r^ 3> 1/to S> 1/tq, so that |^4| S> We have to 
keep in mind that A is negative in this limit. Accordingly, we obtain the following approximation: 



Vt[VA + Bi] = X j-\A\ ( -1 + \/l + B 2 /A 2 ) ~ ^=~^. (A6) 



2' 'V v 'J 2^\A\ ri 



B To 



B. Saha versus Rate Equations 

Some aspects of the baryon-electron system dealt with the recombination rate have been 
discussed in the previous papers (Paper I, II). In Paper I, the perturbed rate equation for hydrogen 
has been presented, while the helium fraction has been neglected and the time-evolution of the 
equation has not been solved. In the present paper we derive a new formula for the perturbation 
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of the fractional ionization taking the helium fraction into account. The evolution of the fractional 
ionization x e is described as (Peebles, 1968) 



d _ {-, Vp 



d -x e = r e n b \l- — )\x e 



2 



2 



1 - y P - (i - y P /2)x s 



C, (Bl) 



where r e is the recombination coefficient, x s is the fractional ionization derived in use of the Saha 
equation and C is the suppression factor. We employ the fitting formula for the recombination 
coefficient by Pequignot et al. (1991) as 

r e = 10 -l3_^? cm 3 s -l (B2) 

where the fitting constants are a = 4.309, b = —0.6166, c = 0.6703, and d = 0.5300, and T4 is the 
baryon temperature in the unit of 10 4 -ff, i.e., T4 = T\ } /10 4 K. From the Saha equation, x s can be 
derived as 



(1 - y p /2fx 2 s = (27rm e T)§ c _ 136eV/T ? ^ 



3 

) 

1 - Vp ~ (1 - y P /2)x s ra b (2vr) 3 " 
The suppression factor C in the right hand side is 



1 + ^ An is 

l + K(A + /3 e )n ls ' lC4J 

where ni s is the number density of hydrogen in the electron state, A is the two-photon decay rate 
from the excited state to the ground state, f3 e = r e (m e T/27r) 3 / 2 e _3 ' 4eV/ ' T and K = (a/a)A 3 1 /87r 
with X a being the Lyman alpha photon wave length. 

Let us derive the perturbation equations on very small scales. We introduce the perturbed 
variables, 

x e = x e {rj) + 5x e (r],x), (B5) 

x s = x s (rj) + 6x s (rj,x), (B6) 

nh = n b (r?)(l + <5 b (r?,x)), (B7) 

where ~ denotes the background quantities. We ignore the temperature perturbation 5T since we 
are only interested in perturbations on very small scales where ST is erased by Silk damping. 



The unperturbed quantities x e {rj) and x s (rj) satisfy the same equations as (plf ) and ( |B3|) with 
replaced rib, x e , and x s by x e , and x s , respectively. Then, from the Saha equation (B3), we 
find the linear perturbation equation for 5x s as, 

2(1 - y p ) - (1 - y p /2)x s K ' 

No perturbation of the fractional ionization can be generated before recombination since 
x s = (I- 2/ P )/(l - 2/ P /2). 



From equation (Bl), on the other hand, we can obtain the linear perturbation equation for 

oXg as, 



- Jt 5x * 



Cr e n h (1 - y p /2) < 5 h x 2 e + 5x e 



2x e + 



(1 " Vvl^l 



1 - Vp - (i - y P / 2 ) x s 



5C 
+ IT 



[1 - y P - (1- 2/ P /2)a; e ] — 



where 
5C 



2/p - (1 - y P /^)xs 
K(3 e n h [(l - y p - (1 - y p /2)x e )5 h - (1 - y p /2)fe e ] 



C7 [1 + £TAn b (l - y p - (1 - y P /2)x e )][l + if (A + /3 e )n b (l - y p - (1 - y p /2)x e )} ' 



(B9) 



(BIO) 



Fig. 7(a) plots the fractional ionization x s and x e for the fiducial CDM model (Oo = 1.0 and 
h = 0.5) and the A CDM model (Q,q = 0.3 and h = 0.7) with high and low baryon densities, i.e., 
n h h 2 = 0.01 and 0.02. 

The evolution of the Fourier coefficients of 5x s and 5x e are plotted in Fig. 7(b). Here we 
adopt the COBE normalized baryon density perturbation with the wave number k = 7120Mpc _1 
which corresponds to 1M & . The nonlinear growth is taken into account to calculate the baryon 
density perturbation. Before recombination, x e ~ x s and 5x e ~ 5x s . As the recombination process 
proceeds, x s rapidly decreases and 5x s falls as rapid as x s because 5x s is roughly proportional to x s 
in equation ( |B8[ ) if x s <C 1. On the other hand, if we employ the rate equation, the recombination 
process proceeds much slower and eventually there remains the residual ionization as is shown in 
Fig. 7(a). This residual ionization results in the residual perturbation of the fractional ionization 
too. Therefore thermal equilibrium is not a good approximation for neither x nor 5x. 

As is discussed in §2, the drag force due to the perturbations of the fractional ionization is 
parameterized by ck/td- Here we define a s = 8x s /{x s 5x>) and a e = 5x e /(x e 5b) to distinguish a's 
obtained by the Saha equation and the rate equation, respectively. 



One can easily derive a s from equation (p38|) as 



(1 - Vp) ~ (1 ~ Vp/tyxs 
2(1 - y p ) - (1 - y p /2)x s 



(Bll) 



Before the recombination epoch, a s = 0. And a s = —1/2 when the fractional ionization can 
be neglected. This value —1/2 is assumed in Shaviv's calculations. On the other hand, a e is 
numerically calculated. Fig. 8 plots a e and a s as a function of (1 + z) for various cosmological 
models. It is found that a e 's are nearly twice as large as q s 's during recombination when the 
nonlinear growth mechanism works. 

In Fig. 9, we plot |a|/rD for same cosmological models. Since the radiation drag force is 
proportional to |q;|/td, we expect larger amplification of the baryon density perturbations if we 
employ the rate equation instead of the Saha equation which Shaviv has used. Moreover, the 
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Fig. 7. — Fractional ionization (a) and its perturbation (b) from the Saha equation and the rate 
equation. We have two values of Qbh 2 ; 0.01 and 0.02 denoted by thin and thick lines, respectively. 
The solid and dashed lines describe the models with = 1.0 and h = 0.5 and = 0.3 and 
h = 0.7, respectively, for the rate equation. In panel (b), long dashed-short dashed lines denote the 
models with = 0.3 and h = 0.7 and dot-dashed lines denote f^o = 1-0 and h = 0.5 for the Saha 
equation. Since the fractional ionization depends only on the value of Q^h 2 in the Saha equation 
(B3), one can find only two lines for x with the Saha equation in panel (a). The wave number k is 
adopted as 7120MPC" 1 . 
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Fig. 8. — Parameter |a| (eq. from the rate equation and the Saha equation for various 

cosmological models. When the fractional ionization is small enough, |a| = 0.5 in use of the Saha 
equation. The wave number k is adopted as 7120Mpc . 
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Fig. 9. — Ratio of \a\ to td from the rate equation (thick lines) and the Saha equation (thin lines) 
for the fiducial CDM models O = 1.0, ft b = 0.04, and h = 0.5 (solid lines) and Q = 1.0, ^ b = 0.08, 
and h = 0.5 (dashed lines). The ratio depends only on Q^h 2 in use of the Saha equation (see eq. 
| Bill and eq. ]17|| ) but also depends weakly on Qq in use of the rate equation. The quasi- nonlinear 
radiation drag force is proportional to |a|/ro. 
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nonlinear mechanism works until later epoch [z ~ 800) with the rate equation due to the residual 
x e and 5x e . 

REFERENCES 

Fukugita, M. k Kawasaki, M. 1994, MNRAS, 269, 563 

Gnedin, N. Y., k Ostriker, J. P. 1997, ApJ, 486, 581 

Jones, B.J.T k Wyse, R.F.G., 1985, Astron. Astrophys. 149, 144 

Nishi, R., Susa, H., Uehara, H., Yamada, M., k Omukai, K., 1998, Prog. Theor. Phys., 100, 881 
Ostriker, J. P. k Gnedin, N. Y., 1996, ApJ, 472, L63 

Peebles, P.J.E., 1993, in Principles of Physical Cosmology (Princeton University Press, Princeton, 
New Jersey) 

Peebles, P.J.E., 1968, ApJ, 153, 1 

Pequignot, D., Petitjean, P. k Boisson, C, 1991 A&A 251, 690 

Sato, H. 1971, Prog. Theor. Phys. 45, 370. 

Shaviv, N. 1998 MNRAS, 297, 1245 

Haiman, Z., k Loeb. A., 1997 ApJ, 483, 21 

Silk, J. 1968, ApJ, 151, 459. 

Weinberg, S. 1971, ApJ, 168, 175. 

Yamamoto, K., Sugiyama, N., k Sato, H., 1997, Phys. Rev. D56, 7566 (Paper I) 
Yamamoto, K., Sugiyama, N., k Sato, H. 1998, ApJ, 501, 442 (Paper II) 



This preprint was prepared with the AAS IATgX macros v4.0. 



